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Abstract. 

A complete solution to the problem of setting up Wigner distribution for iV-level 
quantum systems is presented. The scheme makes use of some of the ideas introduced 
by Dirac in the course of defining functions of noncommuting observables and works 
uniformly for all N. Further, the construction developed here has the virtue of being 
essentially input-free in that it merely requires finding a square root of a certain A^ x A^ 
complex symmetric matrix, a task which, as is shown, can always be accomplished 
analytically. As an illustration, the case of a single qubit is considered in some detail 
and it is shown that one recovers the result of Feynman and Wootters for this case 
without recourse to any auxiliary constructs. 
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1. Introduction 

There has been considerable interest for some time in extending the method of Wigner 
distributions to describe states of quantum systems, originally developed for the case of 
continuous Cartesian coordinates and momenta P-jHl, to the case of finite-dimensional 
quantum systems |I]-[TT]. 

In its original version, the Wigner distribution is a function on the classical phase 
space, real but not necessarily pointwise nonnegative, which describes completely any 
pure or mixed quantum state. Even though it cannot be interpreted as a probability 
distribution on phase space, it does lead to the correct marginal position and momentum 
probability distributions as determined by quantum mechanics. 

Among the early efforts to set up Wigner distribution for states of quantum systems 
with a finite-dimensional state space, one may mention the work of Feynman and of 
Wootters |3]. In the former, attention was devoted to the two-dimensional case, drawing 
on the treatment of spin in quantum mechanics. In the latter it was shown that one has 
to treat separately the cases where the dimension of the state space is a power of two, 
and those where it is odd. In the odd case one has to handle first the case of odd prime 
dimension, and then pass to the general odd situation by forming a Cartesian product 
of the prime cases. 

The approach of Jagannathan 4j on the other hand is based on the Weyl-ordered 
unitary operators for translations on a phase-space lattice, leading to the discrete 
Wigner distribution through the associated characteristic function. The more recent 
independent work of Luis and Pefina ^U] uses a similar approach, but presents a 
thorough analysis of the problem. 

In a previous paper ^2] it has been shown that one can arrive at the Wigner 
distribution concept, in the case of continuous variables, by a novel route starting from 
an idea of Dirac fT^ to describe a general quantum-mechanical operator by a collection 
of mixed matrix elements, using vectors chosen from two different orthonormal bases 
in Hilbert space. The steps that lead from Dirac's starting point to the expression 
for the Wigner distribution, indeed even the introduction of classical phase space 
ideas to describe operators, are particularly transparent and elementary, and they 
automatically ensure the property of correctly reproducing the quantum-mechanical 
probability distributions as marginals. 

The purpose of the present paper is to show that the same approach based 
on Dirac's method can be used in the finite-dimensional well to set up the 

Wigner distribution formalism, incorporating all the desirable features including the 
reproduction of the marginals. It is worth particularly emphasizing that, denoting the 
dimensionality of the state space by N, the present approach works uniformly for all N; 
there is no need to treat separately the cases of a power of two, N an odd prime and 
N an odd number. In the case N = 2, the earlier results of Feynman and Wootters are 
immediately recovered, without having to call upon the specific properties of spin-1/2 
systems. 
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The construction presented here assumes particular significance in the hght of the 
intimate connection between Wigner distributions and mutually unbiased bases J3] as 
was brought out by Wootters and coworkers in a series of insightful papers jTHI • Mutually 
unbiased bases, in turn, are known to be related to questions pertaining to affine planes 
in finite geometries, mutually orthogonal arrays, complex polytopes and finite designs 
[TE] and one expects that the work presented here would provide a new mathematical 
perspective to some of these questions and their interrelations. 

A brief summary of the present work is as follows. In Section 2 we discuss the 
kinematics of N-level quantum systems. In particular we examine the trace of product 
of two operators A and B and show how this can be expressed as a phase space sum 
of products of mixed matrix elements of the operators involved such that it manifestly 
refiects the symmetry under interchange of A and B . This entails introducing a kernel 
whose properties are investigated in Section 3. In Section 4, we show how, by taking 
the 'square root' of this kernel in a certain fashion, one is directly led to the concept of 
a Wigner distribution associated with operators on a A^-dimensional Hilbert space for 
any N. In Section 5, by way of illustration, we consider the case of a single qubit and 
recover known results with economy. Section 6 contains concluding remarks and further 
outlook. 



2. Kinematics of an A^-level quantum system 

We consider a quantum system possessing independent states, so that its state space 
is a complex (finite-dimensional) Hilbert space Ti.^^^ of (complex) dimension A^. We 
select a particular orthonormal basis for written as with g = 0, 1, A^ — 1, 

to be called the set of "position eigenstates" of the system. Then: 

= ^q,q', = 0, 1, A^ - 1; 

N-l /^x 
<?=0 

A general vector E Ti^^^ is described in this basis by a corresponding wave function 
which is an A^-component complex column vector: 

Af-1 N-l 

By means of an A^-point Fourier series transformation we arrive at a complementary 
orthonormal basis of "momentum eigenstates" \p) withp = 0, 1, A^ — 1. The principal 
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equations are: 

7V-1 



N n 

{p\p') = Sp,p', p,p' = 0,h-,N -1, 

N-l (3) 
p=0 

(qlp) = J—e^^'i^l^. 
^ ' ' ^fN 

Now consider a general operator A on liS^^ . Using either the basis {Iq')} or the 
basis for 1-L^^\ it can be completely described by the corresponding N x N square 
complex matrices {q'\A\q) or {p'\A\p). Following the method of Dirac, however, we can 
equally well describe A completely by the collection of "mixed matrix elements" 
we call this a,n N x N "array" rather than a matrix since operator multiphcation is not 
simply the multiphcation of these arrays thought of as matrices. We also notice that 
with the introduction of such arrays the step to a "phase-space" description of A has 
been taken. More precisely, we define the (left) phase-space representative of A by: 

Ai{q,p) = {q\A\p){p\q) = Tr\^A\p){p\q){q\^ = -^{q\A\p) exp {-2mpq/N} . (4) 

(By interchanging the roles of q and p we can equally well define an expression 
Ar{q.p) = {p\A\q) {q\p) , however we will work with the quantities Ai[q,p)). The 
following are immediate consequences of this definition: 

{q,p) = {q\A\q), 

p 

J2A{q:P) = {p\A\p), 
q 

J^Ai (g,p) = Tr{l}. 

<i,p 

We may notice at this point that even for hermitian A, Ai(q,p) is in general 
complex. 

Now take two operators A and B and the trace of their product. We can express 
this in terms of their (left) phase-space representatives as follows: 

Tr[AB] = Nj2MQ,P)Br{Q,P) = Y^Y.(^\A\p){p\q){q'\B\p'){p'\q) 

q,p q,p q',p' 

= EE^' ^) P'^ P'^ (^'' P'^ ' (6) 

q,p q',p' 

where: 

Ki{q,p-<^,p') = N\q\p)ij)\<^)W){p'\q) = exp{27:i{q - (^){p-p') IN} . (7) 
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Thus an important phase-space kernel Ki has been introduced. We note in passing that 
(apart from the A^^ factor) it is a four- vertex Bargmann invariant, so its phase is an 
instance of the kinematic geometric phase [T7] . 

The study of the detailed properties of Ki will lead us to the solution of setting up 
a physically reasonable Wigner distribution, for any value of the dimension N. 

3. Properties of the Kernel Ki 

We can regard Ki {q,p; q' ,p') as defined in Eq. ((Zj) as constituting a complex square 
matrix of dimension A^^, with the first pair of arguments {q,p) being row index and 
the second pair {q',p') column index§. We denote by IC^^ ) a complex linear space of 
dimension iV^, made up of vectors / with components f{q,p)'. 

f e /C(^') ^/(g,p),g,p = 0,1,2 ■■■,iV-l. (8) 
It is to be understood that these vectors are "periodic" in the sense that 

f{q + nN,p + n'N) = f{q,p),n, n' = 0, ±1, ± 2, • • • . (9) 
The norm is defined in the natural way by 

N-l 

ii/ip = (/,/)= E (10) 

q,p=o,i,--- 

Then Ki acts on such vectors according to 

iKif)iq,p) = J2Ki{q,p;q',p') (11) 

q',p' 

The following properties are immediately evident: 

• Symmetry: 

Ki{q,p;q',p') = Ki{q',p';q,p); (12) 

• Essential unitarity: 

J2Ki {q,p; q',p') Ki {q",p"- q',p'Y = iV^ ^ V'S (13) 

q',p' 

• Translation invariance: 

Ki{q + qo,p + po;q' + qo,p' + po) = Ki {q,p; q' ,p') , (14) 

qo,Po = 0,1,2, N -1. 

Here and in the following we interpret translated arguments q + qo,p + Pq, ■ ■ ■ as always 
taken modulo A^, so that they always lie in the range 0,1, ■ ■ ■ , N — 1. Property (|T!^ 

§ We introduce below a more compact efficient notation to express this. 
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means that any eigenvalue of Ki is of the form A^e*"^ for some phase <f. In addition to 
the above, the following 'marginals' properties are also evident from the definition (|7j): 

y^.Ki {q,p] q',p') = NSgqi, independent of p, (15a) 

p' 

y^^Kj {<l,P'i Q ,p) = NSppf, independent of q. {15b) 

q' 

These are particularly important for the Wigner distribution problem, so we explore 
them in some detail and relate them to the eigenvalue and eigenvector properties of Ki. 
From either one of Eqs. we get the (weaker) relations: 

"^^Ki {q,p; q',p') = N, independent of q,p. (16) 

q'p' 

Let us introduce a single symbol a to denote the pair (g, p) by the definition: 

a = qN + P + 1. (17) 

Thus a runs from 1 to A^^: for q = 0,p = 0,1,.. .,N — 1 we have a = 1,2,.. .,N; for 
q = l,p = 0,1, N — 1 we have a = + 1, + 2, 2N; and so on. For summations 
and Kronecker symbols we have the rules: 

9 , siqma=p+l,N+p+l,2N+p+l,...,{N-l)N+p+l 
p fixed 

qp (7=1 
6 era' ^qq' ^pp' • 

We hereafter use a or q,p interchangeably as convenient. The kernel Ki {q,p; q',p') can 
now be written as Ki {a; a'), while vectors / G /C^^ ) have components /(ex). In addition 
to the properties (fT^ . (fT^ . ()15all56|) . (fT^ we have the trace property following from 

TrKi = Y,Ki{a,a) = N\ (19) 

cr 

With this notation one can now see that the marginals properties fll5all56|) can 
be expressed as follows. For each q' = 0, 1, A^ — 1 we define a vector Uq/ in /c(^ ), 
forming altogether a set of A^ real orthonormal vectors (not a basis!) by: 

Ua' (cr) = —;=6aa', independent of p, 

WiV (20) 

{Uq', Uq) = 6q>q. 

Then Eq. pSaj) translates exactly into the statement: 

KiUq = NUq, q = 0,1,..., N-1. (21) 
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Similarly, for each p' = 0, 1, — 1 we define a vector Vp' in /c(^ ) forming altogether 
a set of real orthonormal vectors (again, not a basis!) by: 

Vd' (o") = —;=Svv', independent of o, 
^ y/N (22) 

{Vp', Vp) = 6p'p. 
Then Eq. ()156|) translates into: 

KiVp = NVp, p = 0,1,..., N-1. (23) 
These real eigenvectors Uq and Vp are mutually nonorthogonal: 

{Vp,Uq)=J2Vpi^')UA^') = ^- (24) 

cr' 

This leads to the single linear dependence relation among the 2N (real) vectors Ug, Vp-. 

q p 

which can also be read off from Eqs. ()20|) and ()22j). Therefore, the f/g's and Ip's together 
span an (2A^ — l)-dimensional subspace /C^^^"^) in /c(^ ), over which Ki reduces to N 
times the identity. We can construct an orthonormal basis of {2N — 1) real vectors for 
jq{2N-i) instance by the following recipe: 



q ^ P 

U, = / {Uo + f/i + ... + f/,-1 - jU,) , J = 1, 2, iV - 1, 
V J (j + 1) 

^ 1 (26) 

V, = {Vo + V + ... + V,.i - jV,) ,j = l, 2, ...,N- 1, 

yj (j + 1)^ 

{Uj,U,) = {V,,,V,)=6,'^, (^o,^o) = l, 
{U^,,^,) = {V,,,^o) = {U,,,V,) = 0. 

If the orthogonal complement of /C^^^"^) in /c(^ ), of dimension (A^ — 1)^, is written 
as ^(^-1)', i.e.: 

ici^') = IC(^N-i) ^ j^iN^if ^ (27) 

then we can supplement the basis ()26p for /C^^a^-i) ]^-y (^any) additional real orthonormal 
vectors to span /C^^"^) . The essential unitarity of Ki means that it leaves IC'^^~^^ also 
invariant; the transition from the original (standard) basis of /c(^ ) to the present one 
can be accomplished by an element of the real orthogonal rotation group SO (A^^), thus 
preserving the symmetry ()12j) of Ki. Therefore the matrix Ki has the following structure 
in a (real) basis adapted to the decomposition (j2Z|): 

' A + iB ' ^ 
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The unit matrix is of dimension (2A^ "1); while the two real (A^ — 1) ^-dimensional 
matrices A and B obey: 

= A, = B, AB = BA, 

A' + B' = N'-I^^_,^.^^^^,^2, (29) 

Tr{A} = -N{N -1), Tr{B} = 0. 

Thus the matrix A + iB can definitely be diagonalized by a real rotation in (A^ " 1)^ 
dimensions, i.e., by an element of 5*0 ((A^ — 1)^), and each eigenvalue of A + zB is of 
the form A^e**^ for some angle (p. 

It now turns out that we can carry through this diagonalisation process explicitly. 
The translation invariance ()14j) of Ki means that the eigenvectors of Ki can be 
constructed as "plane waves" in phase space. We can obtain a complete real orthonormal 
set of vectors of Ki in /C^^ ) by this route, recovering the subset of eigenvectors as 
part of a complete set. 

For each point ctq = {qo,Po) we define a unit vector Xao with components 

Xao (^) = ^ (2^^ {qoP + Poq)lN) (30) 

(we see that condition is indeed obeyed). Thus we have exactly A^^ vectors Xo-o- 
Using the modulo A^ rule for phase space arguments we then easily obtain the following: 

Kix., = Are-2-^««P°/^x.o, (31«) 

(Xa^,Xao) = ao- (316) 

Therefore we have achieved full diagonalisation of Ki, with {Xcrg} forming an 
orthonormal basis in /C^^ ). The previously found (real) basis for the subspace IC^'^^~^\ 
made up exclusively of eigenvectors of Ki with eigenvalues A^, is essentially the subset 
of (2A^ — 1) vectors Xqo,o for go = 0, 1, ■ ■ ■ , A^ — 1 and Xo,po for po = 1, ■ ■ ■ , A^ — 1. Indeed 
we find 

N-l 

(32) 

^ N-l ^ ' 

90=0 

The remaining (A^ — 1)^ eigenvectors Xo-o for go, = 1, 2, ■ ■ ■ , A^ — 1 span the orthogonal 
subspace /c(^~^ ). Here we have in detail the following structure. The two eigenvectors 
Xgo,po XN-qo,N-po are degenerate, and their components are related by complex 
conjugation : 

KiXN-qo,N-po = Are-2-«oPo/iv^^_^^^^_^^. (33) 

XN-qo,N-po{(^) = Xqo,PQ{(^)*- 

Therefore we have a pattern that depends on the parity of A^. For odd A^, we 
have (A^ — 1)^/2 distinct degenerate pairs of mutually complex conjugate orthogonal 
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eigenvectors {x<?o,po> XN-go,N-po} for go = 1, 2, ■ ■ ■ , (A^ - l)/2 and po = 1, 2, ■ ■ ■ , - 1. 
For even we have one real eigenvector Xn/2,n/2 with eigenvalue A^(— 1)^/^ followed 
by ((A^ — 1)^ — l)/2 distinct degenerate pairs {xqo,po, XN-qo,N-po} where we omit 
% = Po = N/2. In either case it is clear that by passing to the real and imaginary 
parts of Xgo,po' while leaving Xn/2,n/2 unchanged, we get a real orthonormal basis for 
/(^(^-i) in which the matrix A + iB of Eq. (j^Hj) is diagonal. 

Equipped with these important properties of Ki we turn to Eq. (jH)) from where we 
can find the route to the Wigner distribution. 

4. The Kernel ^ and the Wigner Distribution 

Motivated by the structure © for Tr fo^ ^"wo general operators A and B on 

we try to express the kernel Ki {q,p; q',p') in the form: 

Ki{a,a') = J2a^",'^)acr",a'), (34) 

cr" 

with suitable conditions imposed on ^. The desirable conditions are, as with Ki itself: 
symmetry, essential unitarity, translation invariance and marginal conditions similar to 
Eqs. (Il,^all5/)|l for Kr. 

e(a,a')=e(^x',a), (35a) 
^e(a,a')e(^",0* = A^5.a", (356) 

cr' 

^{q + qo,p + Po;q' + qo,p' + Po) = ^{q,p;q',p'), (35c) 

^U, = VNUg, ^Vp = VNVp. {35 d) 

Here we have expressed the last marginals conditions already in terms of the eigenvectors 
(not all independent!) Ug, Vp of Ki lying in }C^'^^~^\ More explicitly they read 

X^^(?,p; q',p) = "/NSqq., 
<?' 

The detailed analysis of the eigenvectors and eigenvalues of Ki in the previous section 
immediately leads to solutions for ^. The translation invariance of ()35c|) is ensured by 
arranging that the "plane waves" eigenvectors Xgo,po eigenvectors of ^ as well. 

We take ^ to obey: 

^ Xgo,o = v^Xgo.o, go = 0, 1,---,A^- 1; (37a) 

^ Xo,po = ^Xo,po, Po = 0,1,---,A^-1; (37&) 

i (XgcPO or XN-qo,N-po) = ± ViV e-'^^^P"/^ (Xgo,po XN-qo,N-po) , (37c) 

qo,Po = l,---,N -1. 
In the subspaces /C(2A^-i) q^^^ /c(^~^ ) we then have: 
e = ViV-I on/C(2^-i), 

e = (A + zB)^/^ on ^(^-1)1 ^^^^ 
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Equations (|37aj) ensure the validity of the marginals properties ()354 or while Eq. 
()35?)|) is obeyed by construction. It is the symmetry requirement ()35fl|) that dictates 
that in the case of degenerate orthonormal pairs of Ki eigenvectors {xqo,po, XN-qo,N-po} 
we choose the square root of the eigenvalue Ne~'^'^'^'^°^°/^ of Ki in the same way; this is 
expressed in Eq. ()37c|l . Thus we see: for odd N there is a 

2(Ar-i)V2_ fold freedom in the 

choice of ^] for even there is a 

2({JV-i)2+i)/2_foij freedom. In each particular 
square root of A + iB is involved in (jHHj) . 

With any such ^, we can return to Eq. (jHl) and write it in a manifestly kernel- 
independent manner: 

Tr [ab] = nJ2^ {q,p) B iq,p) , (39) 



where: 

A (g, p) = (g, p; q', P) Ai (g', p') (40) 



q',p' 



^X]^ W\A\p'){p'\q'), 



N 

q ,p- 

with a similar expression for B{q,p) in terms of B. We will show below that for 
hermitian A, A{q,p) is real. Combining Eqs.© and ()35dj) we have ensured the marginals 
properties: 

Y,A{q,p) = {q\A\q), 

A - (41) 

{q,p) = {p\A\p). 

q 

For the density matrix p describing some pure or mixed state of the iV-level system, 
we then have the real Wigner distribution: 

W{q,p) = -l=Y^^{q^p-q\p') {q'\p\p'){p'\q') (42) 



q',p' 

and by Eqs. 1)411) the two marginal probability distributions in q and p are immediately 
recovered. In particular, we find that for position eigenstates and momentum eigenstates 
the freedom in the choice of ^ ( which in any case is limited to its action on /c(^~^ )) 
does not matter and we get the anticipated results: 

p = \q'){q'\ =^ W{q,p) = — 6qq/, independent of p, 

^ (43) 
p= \p'){p'\ =^ W {q,p) = —Sppf, independent of q. 

Returning to the Wigner distribution (|42p we may rewrite it as 

W{q,p) = j^Tr{pW{q,p)} (44) 

by introducing elements of Wigner basis jlHI or phase point operators [12] : 

W{q,p) = v^5^e(g,P;g',p') {p'W)\p'){q'l (45) 

q',p' 
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It is an interesting exercise to verify, by combining the definition (|7j) of Ki and (j34l35ftp , 
that these are hermitian: 

W{q,py = W{q,p). (46) 

This proves that W{q,p), and more generally A{q,p) = Tr{AW{q,p)} for hermitian A, 
are both real. In addition one can check, by virtue of eqs ()34p35al -(i). they satisfy 

Tr{W{a)} = 1, 

Tr{W{a)W{a')} = N5,^>. 



(47) 



5. The case of = 2: the Qubit 



This case is particularly interesting in that earlier treatments have had to treat it on 
its own, in a sense in an ad hoc manner, as distinct from an odd prime or an odd 
integer. In the standard basis for the two-dimensional Hilbert space H^'^^ made up of 
|g) for g = 0, 1, accompanied by its complementary basis \p), the matrix Ki {q,p] q',p') 
is the following: 

-1 \ 
1 



-1 



-1 



\ 



-1 



(4^ 



The rows and columns are labelled in the sequence: (g,p) = (0, 0) , (0, 1) , (1, 0) , (1, 1), 
and the matrix elements are read off from Eq. (jZj). The three orthonormal eigenvectors 
of Ki with eigenvalue 2, spanning the subspace /C^'^^ of the general treatment in Section 
3 are: 

/ 1 \ 



" 2 



/ 1 \ 

1 
1 

VI/ 



Ui = - 
2 



-1 



2 



/ 1 \ 

-1 



1 



(49) 



We choose the fourth eigenvector of Ki, 
TrKi = 4, to be: 

/ 1 \ 



^-2 



-1 
-1 

V 1 / 



with eigenvalue necessarily —2 since 



(50) 



Then the kernel ^ can be immediately synthesized from: 

and in the standard basis turns out to be: 

f 3 + i 1-i 1-i -l + i\ 

1 1-i 3+i -1+i 1-i 

2V2 1-i -1 + i S + i 1-i 

\ -1 + i 1-i 1-i 3 + i J 



(51) 



(52) 
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Using the above matrix elements of ^ in (j45|) we obtain, for the phase-point operators: 








2 " / \ 2 




(53) 



>v(i,o)= I I , w(i,i: ' 



=1^ / ' ^^'^^ I ^ 1 



and thereby recover the results of Feynman and Wootters |3] and hence also the 
connection between sums of phase point operators along striations of the qubit phase 
space [T3] and the mutually unbiased bases for = 2. 

For the density operator p = |(I2 + a • a), a ■ a < 1 describing a general state of a 
qubit one can easily calculate the corresponding Wigner distribution using or (jUj). 
The results arranged in the form of a matrix read: 



1 + 01 + 02 + 03 1 + 01—02 
1 — Oi — O2 + O3 1 — Oi + O2 



(54) 



Using this result it is instructive to verify the validity of for A = pi = |(I2 + a ■ a), 
and B = p2 = |(I2 + b ■ cr). Further, it is easily seen from (j5H) that the maximum 
positive and maximum negative values of the Wigner distribution for a qubit occur 
when |oi| = I02I = I03I = l/^/S. 

As a final remark, we notice that in calculating the kernel ^, square root of Ki, 
we might have chosen W in Eq. ()5H) to be the eigenvector corresponding to the 
eigenvalue —i\pl instead of +i\/2. This results in changing +i{—i) with — i(+i) in Eqs. 
and thus in an interchange of the role of W (0, 0) and W (0, 1) with W (1, 0) and 



W (1, 1) respectively. Correspondingly, the coefficient 02 in Eq. would change sign 
everywhere. This however can be of no physical consequence as is reflected in the fact 
that the marginals obtained by summing over g or p are independent of 02. 



6. Concluding remarks. 

To conclude, we have developed a method of constructing Wigner distribution for 
A^-level systems which is remarkable in its directness and economy and works uniformly 
for all A^. The construction is entirely algebraic and solely involves finding a square root 
of a certain A^^ x A^^ complex symmetric matrix. No other auxiliary inputs are required. 
As an illustration, we have worked out the qubit case in some detail and obtained results 
already known in the literature in an extremely economic fashion. The construction 
presented here provides a fresh perspective to several questions pertaining to quantum 
tomography in finite state systems and to those associated with finite geometries which 
are currently being investigated with great vigour owing to their relevance to quantum 
information theory. We hope to return to some of these elsewhere. 
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